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We present solutions of the time dependent Schrfidinger equation for a SQUID ring coupled to 
an electromagnetic field, both treated quantum mechanically. We that show the SQUID ring can 
be used to create a maximally entangled state with the em field that without dissipation remains 
constant in time. Using methods familiar in quantum optics, we extend the model to include the 
effects of coupling this system to a dissipative environment. With this model we show that although 
such an environment makes a noticeable difference to the time evolution of the system, it need 
not destroy the the entanglement of this coupled system over time scales required for quantum 
technologies. 
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Experiments to probe the quantum properties of su- 
perconducting Josephson weak link (SQUID) circuits [ij, 
S [10 (in this work a superconducting ring enclosing 
a single weak link) have highlighted the possible use of 
such devices in emerging quantum technologies such as 
quantum computation, communication and information 
processing In an earlier paper, motivated in part 

by experiment we solved the time dependent Schrodinger 
equation (TDSE) for an isolated SQUID ring modulated 
by a classical, monochromatic, microwave field [fj. Re- 
cent studies of a fully quantum mechanical model 
of electromagnetic (em) fields coupled to SQUID rings 
have, from a theoretical standpoint, furthered our un- 
derstanding of the manner in which such circuit may be 
utilised. 

An essential feature revealed by these studies is 
that a quantum mechanical SQUID ring is highly non- 
perturbative in nature, a property that can be used to 
affect very strongly the non-linear dynamical behaviour 
of SQUID ring based quantum circuit systems. Further- 
more, these effects can be controlled through the external 
(bias) magnetic flux $ x applied to the ring. Even so, al- 
though several interesting new results came from the in- 
vestigation of coupled ring-em held systems, in each case 
we chose, for ease of computation, a SQUID ring Hamil- 
tonian description which was time independent. How- 
ever, there are many circumstances where time depen- 
dence needs to be introduced into this description. In 
the work described here this arises because we wish to 
entangle a SQUID ring with an em held mode. As we 
will show, this can be achieved by applying a time de- 
pendent magnetic flux to the SQUID ring which changes 
slowly enough to maintain the ring-held system in its 
current state but fast enough to ensure strong entangle- 
ment. In this paper we solve the TDSE or, more gener- 



ally, the master equation where dissipation involved, for 
tensor products of SQUID ring-em field systems with a 
particular time dependent potential. We show that we 
can use the flux dependent, non-perturbative nature of 
the SQUID ring to great utility in controlling its level of 
entanglement with an em field. As we shall also demon- 
strate, the level of this entanglement can decay with time 
or remain constant depending on whether or not dissipa- 
tion (decoherence) is present. 

In dealing with a SQUID ring coupled to an em field 
as quantum mechanical objects we represent the field as 
an equivalent circuit, a parallel, inductor-capacitor (LC), 
oscillator |l(J . We also assume that the couplin g be tween 
the SQUID ring and the em field is inductive 12]. The 
schematic for this coupled circuit system is shown inset 
in[U 

With reference to the inset in figure ^ we can write 
down a Hamiltonian for the coupled two mode (SQUID 
ring + em field oscillator mode) system of the form 

H = H e + H s — H es (1) 

where H e and H s are, respectively, the Hamiltonians for 
the em field and the SQUID ring and the interaction en- 
ergy between them (for our purposes purely inductive) 
is H es . For these we adopt the convention that oper- 
ators with the subscripts s or e represent, respectively, 
those quantities associated with the SQUID ring or the 
em field. We also note that we use the first four energy 
eigenstates of each Hamiltonian to represent the individ- 
ual circuit component of the system. This number was 
found to be sufficient to compute accurately the results 
presented in this paper. 

The Hamiltonians for these quantum circuits can 
be expressed in terms of the conjugate variables of 
magnetic flux $ (coordinate) and electric displacement 
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flux Q (— > — ihd /d&)&nd the appropriate circuit param- 
eters. Thus, for a SQUID ring (inductance A s , weak 
link capacitance C s ) $ is the total magnetic flux thread- 
ing the ring and Q is the total electric displacement 
flux (displacement charge) between the electrodes of the 
weak link; for the field oscillator these flux and charge 
variables relate, respectively, to the magnetic and elec- 
tric displacement fluxes threading the equivalent induc- 
tor (inductance A e ) and capacitor (capacitance C e ). In 
each case these operators satisfy the usual cannonical 
commutation relation [<f>, Q] — ih. We assume that 
both the SQUID ring and em field oscillators (frequen- 
cies lo s /2tt (= l/2wy/A s C s ) and u> e /2n (= l/2ny/A e C e ) , 
respectively) operate in the quantum regime, i.e. such 
that Huj s , hio e ^ k^T for temperature T. 

Given a SQUID ring inductance of A s and a weak link 
capacitance of C s , the SQUID ring Hamiltonian is, as 
usual, translated in external bias flux and takes the 
time dependent form 
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where hv/2 — 0.0215 &^/A s is the Josephson phase 
coherent coupling energy, $o (= h/2e) is the supercon- 
ducting flux quantum and we note that the last term, 
Qs-§i^x{t) in @ arises from our translating the Hamilto- 
nian by *& x {t) (see for more details). For this quantum 
regime we set C a = 1 x 10~ 16 F and K s = 3 x 10" 10 H, 
parameter values which are appropriate for a mesoscopic 
SQUID ring operating at a few K [Icj . 

The Hamiltonian for the em field mode, which is some- 
what simpler but equivalent to that for a simple harmonic 
oscillator, can be written as 



H e = 
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while the energy associated with the inductive coupling 
between the ring and the em field mode is given by 



H es = ^$ s $e 

A., 



(4) 



where /i es = 0.01 is the ring-field magnetic flux linkage 
factor. For simplicity, and to facilitate the computations, 
we set C e = C s and A e = A s throughout this work. 

Using Q), we showed in our previous work (l(3 | that 
energy could be exchanged between the quantum modes 
of a SQUID ring-em field system around particular values 
of the static flux applied to the ring. For the SQUID ring 
(i = s) and the em field (i = e) the time averaged energy 
expectation values are found from 



((Hi)) = lim - 

r — *oo t 



Examples of exchange (or transition) regions in these 
time averaged energies are shown in figure ^ around the 
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FIG. 1: Time averaged energy expectation values for the ini- 
tial state | l e s ) for (a) the em field and (b) the SQUID ring 
as functions of $1. 

arrowed point A (i.e. where the energy difference be- 
tween two of the ring energies is very close to nhuj e , n 
integer). As in the rest of the paper, the initial state was 
taken to be \ip (0)) = |l e s ), i.e. with the cm field mode 
in its first excited Fock state and the SQUID ring in its 
ground energy state. In computing these averaged en- 
ergies we clearly could not integrate for the system over 
an infinite time; instead we calculated the time averaged 
expectation values for a time r such that ((Hi)) did not 
change noticeably if r was further increased. 

The actual value of $ x around which exchange regions 
develop depends on circuit (ring and field oscillator) pa- 
rameters and the initial condition of the system. For 
example, for the particular case presented in figure ^ the 
point A, at $ x — 0.42864 $ (together with its twin at 
Q x = 0.57136<i>o, as shown), is where significant energy 
is exchanged between the SQUID ring and the em field. 

The energy exchange, which reaches a maximum at the 
centre of the exchange region, follows the strength of the 
coupling between the SQUID ring and the em field mode. 
Other subtle properties of this two component quantum 
circuit system also follow this changing coupling strength, 
in particular the quantum entanglement between the ring 
and the field mode. Entanglement plays a central role in 
considerations of quantum technologies 0, where it 
is often of the first importance that strong entanglement 
is maintained as the quantum system evolves with time. 
With this in mind, and as the central result in this paper, 
we set about trying to create a state of maximal entan- 
glement that remains constant with time. In our chosen 
system this takes the form 

|^(t)) = i(e^|l e s ) + e^|0 e l s )) (5) 
To implement this we now make the external flux time de- 
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FIG. 2: (a) Probability of the system being in state |l e s ) 
or |0 e l s ), i.e. |l e s ) or 9* |0 e l a ) (b) system entanglement 
entropies I s and I e (here equal) evolving at point A for initial 
condition |l e s ). 

pendent with a functional form (see inset in figure El (b)) 
given by 

!A t <t 

^^(t-to) t <t<t +t r (6) 
B t>t + U 

where to = 326WJ 1 and the ramp time is t r — I6.6W7 1 . 
Following this change in bias flux, the operating point 
has moved from A to B (& x = 0.38$o) in figure ^ well 
outside the exchange region. The cut off time to is the 
time at which the probability of the system being in state 
|l e s ) is equal to the probability of it being in state |0 e l s ). 
For our two component system we found these probabil- 
ities by solving the TDSE with the ring and field mode 
parameters given above. As there is no exchange of en- 
ergy between the ring and the field mode at point J5, 
it seemed reasonable to assume that at this flux bias the 
system would maintain equal probability of being in each 
of the states |l e s ) and |0 e l s ) at times after to. The 
results of our calculations of these probabilities are pre- 
sented in figure Eta), plotted against dimensionless time 
uj e t. We see from these computed solutions that our as- 
sumption is essentially correct, at least for the particular 
system example being studied. 

Given this prediction of the behaviour of the circuit 
system, it is important to establish a measure of the level 
of entanglement between its two components. In this 
work we adopt the entropic measure of entanglement first 
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FIG. 3: For comparison with figure|21(a) probabilities of being 
in the state |l e s ) or |0 e l a ) and (b) entanglement entropies 
Is and I e for the system evolving in a dissipative environment 
at point A with the initial condition |l e s ) . 



introduced by Adami and Cerf .13], i.e. 

I i = S(p)-S( Pi ) 

where S (p) = Tr (p In p) and p is the density operator. 
If It is negative this implies that the z th component of 
the system is entangled with the rest of the system, an 
inequality that holds even in the presence of dissipation. 
This computed quantity for our system is shown in fig- 
ure [21 (b) as a function of ui e t. As we might suspect, to 
is also the time when the entanglement between the em 
field and the SQUID ring reaches a maximum. From the 
figure we see that after approximately a time to both the 
ring and field mode probabilities have reached values of 
0.5 (figure EJa)) which then remain constant with time. 
Similarly, in Efb) beyond to the entanglement is strong 
(at 0.7) and also constant with time. 

In the above discussion we have seen from solutions of 
the TDSE that we can use a time dependent flux ©, ap- 
plied to the SQUID ring, to make a maximally entangled 
state of the form (J5J. Nevertheless, for this entanglement 
to be of use technologically it must be maintained over 
some reasonable period of time in the presence of dis- 
sipation. In order to estimate the effects of dissipation 
we adopt the open systems approach usually used in the 
field of quantum optics. In this approach the evolution 
of the system is found by solving a master equation for 
the density operator. The master equation for a system 
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FIG. 4: For comparison with figure|5|(a) probabilities of being 
in state |l e s ) or |0 e la) and (b) entanglement entropies I s and 
I e in a strongly dissipative environment at point A with initial 
condition |l e s ). 

subject to thermal dissipation is given by 

= [H , p] + ^2 (Mi + 1) [2a. lP a\ - a\a iP - pa\a, 

i=e,s 

+ l^Mi {^alp^ - aia\p - penal J 

where Tj is the damping rate of each system to the ther- 
mal bath and the mean photon number Mj = 2.746 x 
10~ 5 in each component is related to the temperature 
Th i = 4.2K and frequency uj^ = uj s of each decohering 
bath via Mi — [expQuo^/kBT) — 1] . 

The effect of dissipation on the probabilities and en- 
tanglement entropies of our chosen system are shown in 
figures [3] and for the cases of weakly and strongly dissi- 
pative environments, respectively. Thus, in figure |3| at a 
7 value of 1 x 10^ 5 w s , the effect on the probablities is neg- 
ligible over a few to while the entanglement can be seen 
to weaken only gradually over the same time period. For 
this degree of decoherence the corrsponding decoherence 
time is of the order of 20ns, which is in line with recent 
experimental results on superconducting weak link cir- 
cuits (albeit at much larger weak link capacitances 
(w 10~ 12 F) and lower temperatures (s» 25mK)). This 
suggests that at this level of dissipation the coherence of 
the coupled system could be used effectively for quantum 
coherent circuit operations. It is also apparent from the 
second example of the effect of a dissipative environment 
that substantially greater values of 7 (= 1 x 10~ 4 w s ) i.e. 



stronger dissipation, still do not entirely suppress entan- 
glement over this same time scale of a few to. As can be 
seen, at this higher level of dissipation our system starts 
to become mixed after a sufficient evolution in time. Con- 
comitantly, a time independent, maximally entangled, 
state is not maintained. Clearly, this result suggests a 
practical limit to the level of dissipation (decoherence) 
that can be tolerated if quantum coherent operations are 
to be performed successfully. With this in mind, we note 
that these results do not necessarily imply that dissipa- 
tive effects will be significant, although when designing 
real circuits environmental effects must be accounted for 
in its design. 

In this paper we have shown that, via a time dependent 
external bias flux, a SQUID ring can be used to control 
the quantum state of a coupled em field mode-SQUID 
ring system, including the degree of entanglement be- 
tween these two circuit components. As a corollary of 
this, we showed that the level of this entanglement could 
be set and maintained, with only a slow degradation in 
time, when the system was coupled to a weakly dissi- 
pative external environment. With the current interest 
in using superconducting weak link circuits for quantum 
technologies, this work provides some guide to the time 
scales over which it may be possible to perform quan- 
tum coherent operations in the presence of decohering 
environments. 

We would like to thank Dr. T.P. Spiller for interesting 
discussions on entanglement in quantum circuits. 
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